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FOREWORD 

This  Interim  Technical  Report  covers  research  initiated  by  the 
Aeronautical  Research  Laboratories,  Office  of  Aerospace  Research,  United 
States  Air  Force,  on  Contract  No.  AF  33(615)-2764,  "Research  in  Numerical 
Analysis.”  TCie  work  was  done  under  the  cognizance  of  the  Applied  Mathe¬ 


matics  Laboratory  with  Mr.  H.  E.  Fettis  acting  as  Project  Director. 


ABSTRACT 


In  this  report,  ve  investigate  the  application  of  Miller's  algorithm 
for  backward  recursion  to  the  difference  equation 
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Cv(n)y(n+--) 


=  0 


c0  »  1  ,  Ca(n)  /  0 


a  >  1  , 


n  -  0,1,2,...,  where  Cv(n)  has  a  representation  ac  an  asymptotic  series 
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Kv/«i 


~  n 


-l/to 

V  +  cl,vn  + 


♦1  ,  n- 


and  Kv  is  an  integer,  u.‘  is  an  integer  ^  1  ,  cQ^v  ^  0  unless 
Cv(n)  E  0  .  Also  we  discuss  the  convergence  of  several  methods  which  are 
related  to  Miller's  algoiithm  and  which  con  be  used  when  Miller's  algorithm 
does  not.  converge. 

A  number  of  examples  are  included. 
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I 

SUMMARY  AND  INTRODUCTION 

Computation  by  the  use  of  difference  equations  in  the  backward  direction 
was  introduced  by  J.C.P.  Miller  (1952).  In  this  reference,  he  applied  the 
method,  now  sometimes  called  Miller's  algorithm,  to  the  calculation  of  Bessel 
functions,  and  it  proceeds  essentially  as  follows.* 

Consider  the  difference  equation 

y(n)  -  y(nU)  -  y(n*2)  =  0,x>0,n£0,  (l.l) 

which  is  satisfied  by  the  modified  Bessel  functions  In(x)  and  (-)nKn(x)  .** 
Let  m  be  an  integer  £  0  .  Put 

Am+1(m)  =  0  '  =  L  >  (l*2) 

and  calculate  An(m)  fo’*  0  ^  n  s  m-1  from  (l.l),  i.e., 

An(m)  =  An+l(m)  +  An+2(m)  ,  0  s  n  s  m-1  .  (1.3) 


*  Miller  has  stated  that  he  first  used  the  method  as 
lion  of  Airy  integrals,  see  Miller  (L91G). 

**  All  special  functions  in  this  work  are  defined  as 


an  aid  in  the  compu ta¬ 
in  Krdclyi  et  si  (1953). 
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is  known  (Erdclyi  et  al  (1953),  v.  II,  p.  7). 

Let 

0/2]  k 

n(m)  =  ]>_  (-)  ek^2k^m^  >  (^*5) 

k-0 

where  (jn/2]  means  the  largest  integer  not  greater  than  m/2  .  Then,  by 
using  the  known  asymptotic  properties  of  In  ,  Kn  for  large  n  ,  one  can  show 
that 

lim  An(m)/n(m)  3  I^x)  ,  n  *  0  ,  x  >  0  .  (3,6) 


In  fact,  the  asymptotic  estimates 

In(x)  =  [W.n"1)]  ,  (-)\(x)  3  (-g^  r(n)  [l+OCn-1)]  (1.7) 

follow  from  the  ascending  series  represents  ions  of  In  and  Kn  ,  and,  since 
An(m)  satisfies  (l.l),  it  can  be  represented  as  a  linear  combination  of  the 
(linearly  independent)  solutions  (1.7),  see  the  Appendix.  This  means  that 

An(m)  »  §1(m)ln(x)  +  52(m)(-)nK^(x)  .  (1.0) 


From  (1.2)  and 


(*)Vi<x>  f  l/x  ' 


(1.9) 


wc  conclude  that 
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§l(m)  »  xK^x)  ,  ?-(m)  »  x(..)mV1(x)  . 


(1.10) 


Thus 


An(m)  =  m.'(2/x)min(x)  [l  +  0(m"l)J  ,  m->®  , 
0(m)  =  m!(2/x)m  j~l  +  0(iu-1)  J  ,  m->«  , 


(1.11) 


£30  (1.6)  follOWS. 

The  above  analysis  shows  clearly  why  the  process  converges,  and  also  why 

n 

it  converges  to  In  and  not  to  (-)*Kn:  In  is  very  snail  compared  to  Kn 
as  n-»®  .  Th^s  characteristic  of  Miller's  algorithm,  namely,  that  the  solu¬ 
tion  of  the  difference  equation  to  which  the  process  converges,  if  it  converges, 
must,  in  a  certain  sense  be  the  smallest  solution,  remains  true  when  the  algo¬ 
rithm  is  applied  to  general  homogeneous  difference  equations,  see  our  Theorem 
4.3. 

A  remarkable  feature  of  the  Miller  algorithm  is  that  no  tabular  values  of 
ln  are  needed  in  the  computations,  only  a  normalization  relationship,  such  ae 
(1.4).  Tabular  values  vould  be  required,  of  course,  if  (l.l)  were  used  in  the 
forward  direction,  and  moreover,  when  (l.l)  is  used  in  the  forward  direction 
to  compute  lu  starting  with  initial  values  of  Ic  and  J^,  those  small  errors 
inevitably  introduced  Ln  the  course  of  the  computation  grow  rapidly  with  r.  . 
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5 Vf. 


Such  a  phenomenon  is  called  instability.* 

The  method  proposed  by  Miller  created  enormous  interest,  and  a  number  of 
papers  subsequently  appeared  in  which  the  writers  either  further  treated  the 
application  of  the  method  to  Bessel  functions,  or  else  showed  that  the  method 
could  be  used  to  compute  other  special  functions.  Stegun  and  Abramowitz 
(1957),  Rondels  and  Reeves  (1950),  Goldstein  and  Thaler  (1959),  Corbato  and 
Uretsky  (.1959),  and  Makinouchi  ( 1965a, b)  all  treated  the  computation  of  Bessel 
functions.  Rotenberg  (19GG)  showed  how  the  algorithm  could  be  used  to  compute 
toroidal  harmonics  (i.e,,  Legendre  functions)  and  Miller  himself  applied  the 
method  to  parabolic  cylinder  functions  (1964) . 

Gautschi  (19Gla)  discussed  the  computation  of  repeated  integrals  of  the 
error  function 


00 

y(n)  a  in  erfbx  =  ^/^JrTnlJ  J*  (t-x)ne_t  dt  ,  n  ^  0  ,  (1.12) 


which  satisfy 


*  The  computation  of  (-)nKn  by  using  (l.l)  in  the  forward  direction  with 
initial  values  of  K0  and  Kj_  is  stable,  i.e.,  I'andom  errors  intro¬ 
duced  during  the  computations  do  not  grow  with  n  .  In  general,  a  differ¬ 
ence  equation  can  be  used  efficiently  in  the  forward  direction  only  to 
compute  the  "largest"  solution  of  the  equation.  However,  the  analysis  of 
the  forward  procedure  is  rather  less  of  a  problem  than  the  analysis  of 
Miller's  algorithm,  see  Gautschi  (ca  1962),  and  will  occupy  none  of  our 
attention  here. 
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y(n)  -  2xy(n+l)  -  2(n+2)y(n+2)  =  0  , 


(1.13} 


and  in  a  later  paper  (1961b)  discussed  the  computation  by  backward  recursion 
of  a  number  of  other  functions  defined  by  definite  integrals. 

The  Miller  algorithm  can  be  applied  to  problems  other  than  the  computa¬ 
tion  of  the  special  functions.  Recently,  it  has  been  employed  in  such  diverse 
problems  as  the  calculation  of  successive  derivutives  of  £f(z)/zj  ,  where  f 
is  an  arbitrary  analytic  function  (Gautcchi  (1966))  u.nd  the  computation  of 
coefficients  for  the  Chebyshev  polynomial  expansions  of  functions  vhich  satisfy 
differential  equations  with  polynomial  coefficients  (Clenshaw  (1957), (1362) ) . 

Of  course,  any  numerical  technique  of  such  general  applicability  demands 
a  thorough  theoretical  investigation.  Gautschi  (1561b),  who  analyzed  its  con¬ 
vergence  when  applied  to  an  arbitrary  second  order  difference  equation,  seems 
to  have  been  the  first  writer  to  discuss  the  Miller  algorithm  from  a  general 
point  of  view.  He  continues  this  analysis  in  two  unpublished  works  (ca  1962, 
1963)  using  as  his  main  tools  the  theory  of  continued  fractions  and  the  classi¬ 
cal  asymptotic  theory  of  linear  difference  equations  (e.g.,  the  theorems  of 
Poincare,  Perron,  and  Kreuser),  and  he  applies  his  findings  to  the  computation 
of  Bessel  functions,  Legendre  functions,  the  incomplete  Beta  function  and  the 
numerical  computation  of  Fourier  coefficients. 

By  now  a  great  deal  is  knovn  about  the  application  of  Miller's  algorithm 
to  second  order  difference  equations.  Conditions  on  the  solutions  of  the  equa¬ 
tion  which  will  guarantee  the  convergence  of  the  algorithm  have  been  given  by 
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Gautschi  (1961b)  and  by  Olver  (1964).  Methods  of  increasing  the  efficiency 
of  the  algorithm  by  using  the  adjoint  equation  to  generate  auxiliary  sequences 
have  been  given  by  Shintuni  (1965).  (Our  Theorem  2.7  is  a  generalization  of 
one  of  this  author’s  results.) 

On  the  other  hand,  little  is  knovm  concerning  application  of  the  algorithm 
to  difference  equations  of  arbitrary  order.  Gautschi  (ca  1962)  has  touched 
briefly  on  the  use  of  difference  equations  of  order  a  >  2  ,  but,  unfortunately, 
the  classical  asymptotic  theory  on  which  his  analysis  is  based  does  not  give 
very  realistic  conditions  for  convei-gence  of  the  algorithm.  More  specifically, 
Gautschi  found  it  necessary  to  assume  the  existence  of  a  fundamental  set  for 
the  equation  (see  the  Appendix)  whose  members  exhibited  radically  dissimilar 
behaviour  as  n-*«  ,  namely  if  jyh(n)j  were  the  set  in  question,  and 

yh(n+i)/yh(n)  ~  V*  h  i  \  /  0  >  va  >  vo-l  >  ...  >  vx  ,  (1.14) 

n-*<"  ,  1  Sh  Sir  , 

then  it  could  be  shown  that  Miller's  algorithm  converged  (to  y-^(n))  ,  provided, 
of  course,  that  a  suitable  normalization  relationship  was  known.  Needless  to 
say,  this  condition  is  excessively  stringent,  at  least  for  a  very  wide  class 
of  difference  equations  (cf  our  Theorem  4.2). 

The  purpose  of  this  work  is  to  examine  the  application  of  Miller's  algo¬ 
rithm,  as  well  as  several  related  algorithms,  to  homogc - 'ous  linear  difference 
equations  of  arbitrary  order  with  coefficients  of  a  fairly  general  type. 
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(The  related  algorithms  are  modifications  of  Miller’s  algorithm  which  can  be 
used  when  the  equation  has  no  smallest  solution.) 

First,  in  Chapter  II,  we  formulate  the  algorithms  as  they  are  applied  to 
difference  equations  with  arbitrary  coefficients,  and  investigate  their  con¬ 
vergence  properties.  Conditions  for  convergence  take  the  form  of  rather 
unwieldy  restrictions  on  the  growth  of  solutions  of  the  equation  and  determi¬ 
nants  involving  them.  To  obtain  more  practical  conditions,  it  is  necessary 
to  restrict  somehow  the  form  of  the  coefficients  in  the  equation.  This  is 
done  in  the  chapters  following  by  requiring  that  the  coefficients  possess 
certain  asymptotic  representations  as  n-*»  .  Ue  are  then  able  to  use  as  our 
principal  investigational  tool  the  analytic  theory  of  singular  difference 
equations  which  was  developed  by  Birkhoff  and  Trjitzinsky.  Virtually  all 
difference  equations  encountered  in  practical  applications  are  of  the  specified 
form,  including  all  equations  with  coefficients  rational  in  n  .  In  particular, 
the  computational  procedures  discussed  by  the  preceding  authors  are  included 
in  our  analysis. 

Chapter  III  is  devoted  exclusively  to  an  asymptotic  analysis  of  the  solu¬ 
tions  of  this  difference  equation,  starting  with  the  above-mentioned  theory  of 
Birkhoff  and  Trjitzinsky.  We  then  prove  two  new  theorems  concerning  the  repre¬ 
sentation  of  those  solutions  whose  growth  can  be  described  with  only  elgebraic 
and  logarithmic  terms. 

In  Chapter  IV,  wc  apply  the  asymptotic  results  of  the  previous  chapter 

to  the  problem  of  determining  more  tractable  conditions  for  the  convergence 
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of  the  algorithms  Given  in  Chapter  II.  One  result  is  that,  for  the  type  of 
difference  equation  considered,  at  least  one  of  the  algorithms  will  converge 
to  a  solution  of  that  equation,  provided  tliat  one  can  find  a  fundamental  sot 
for  the  adjoint  equation  in  vhich  no  more  than  two  solutions  have  the  same 
rate  of  growth  (in  absolute  value)  as  n-* «  .  A  consequence  of  this  theorem 
is  that,  for  a  second  order  difference  equation  of  the  specified  type, 

Miller's  algorithm,  or  a  suitable  modification  of  it,  will  always  converge. 

Chapter  IV  contains  a  number  of  examples,  among  which  are  the  computation 
of  the  integral 


y(n)  =  J  exp  j-ta  -  P(t)j  tndt  ,  n  s  0 


(1.15) 


a  an  integer  >-  2  ,  P(t)  a  polynomial  of  degree  (cr-1)  ,  and  the  computation 
of  a  class  of  hypergeometric  functions. 

In  the  Appendix  of  this  work  are  contained  definitions,  notation,  and 
those  properties  of  difference  equations  which  are  frequently  used  in  the  pre¬ 
ceding  chapters.  (References  to  material  in  the  Appendix  are  preceded  by  an 
A,  as  A. 2,  A-VI,  etc.) 


II. 


In  this  chapter,  we  will  discuss  computation  by  backward  recursion  based 
on  the  linear  homogeneous  difference  equation  of  order  o  s  2 

a 

YL  Cv(n)y(n+v)  -  0  ,  C0  -  1  ,  Ca(n)  f  0  ,  (2.1) 

VsO 

where  n  is  an  integer  S  0  . 

The  first  algorithm  proceeds  as  follows.  Let  m  be  an  integer  i  0  . 

Put 

Vi+a-l(m)  =  ^m+c-2(m)  58  •••  “  ^^(m)  =  ®  >  ^(w)  a  ^  >  (2.2) 

and  calculate  An(m)  for  0  i  n  i  m-1  by  recursion  from 

o 

^  Cv(n) Aj.1.+v(,!')  ”  0  •  (2.3) 

v-0 

Suppose  we  are  given  the  convergent  series  (called  a  normalization  rela¬ 
tionship) 


1 3  Z  viM  (2.*) 

k=*0 

where  y^n)  is  a  solution  of  (2.1), 
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Define 


m 

n(m)  =  Yi  MkM 

k=0 


> 


and 


rn(m)  =  An(m j/n(m)  . 


Definition  2.1 
If 


lim  rn(m)  =  y^n)  ,  n  s  0  , 

m->o 


(2.5) 


(2.6) 


(2.7) 


then  we  say  the  computation  of  y^(n)  by  baclarard  recursion  based  on  (2.1) 
and  (2.4)  converges. 


Let  us  analyze  the  above  algorithm.  First,  by  A- IV  we  note  that  (2.1) 


possesses  a  fundamental  set, 
tion,  we  can  write 


and  since  An(m)  satisfies  the  equa- 


o_ 

An(m)  =  5h(m)yh(n)  *  (2-8) 

h=l 


where  is  independent  of  n  .  By  setting  n  =  m  ,  m  +  1  ,  ...  ,  m  +  a  -  1 
in  (2.8)  it  is  found  that 
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?h(m)  *  a'h(m)/D(m)  / 


(2.9) 


where  ,  D  are  defined  in  A-V,  A-VI.  D  is  not  zero  because 
A-I,  A-V,  and  the  fact  that  Ca(n)  /  0  . 

Thus 


o 

H  Th(m)yh(n) 

h=l _ 

a  m 

Y  Th(m)  21  l*yh(k) 

h=l  k=0 


Theorem  2.1 

Let  T-^(m)  ^  0  for  m  sufficiently  large. 
Define 


rn(n)  = 


and  this  formulation  leads  to 


Rh 


Rh(m) 


V”) 


m 

Sh  s  ShC®)  3  Y  Vh(k) 

k=0 


1  5  h  5 


Nov  suppose 


Urn  Rp,  =  lim  •-  0  ,  2  £  h  £  a  . 
m-»»  m-»® 


Then  the  computation  of  y-j^n)  by  backward  recursion  based  on 
converges. 


of  properties 

(2.10) 


o  .  (2.11) 

(2.12) 

?.l)  and  (2.4) 
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Proof: 


By  (2.8)  and  A-VI  we  have,  for  k  ,  m  >■  0  , 


a 

H  Cff.v(m+v)Ak(»n»)  »  0  ,  (2.17) 

VaQ 


m+cr  a 

H  H  *  0  .  (2.18) 

k=0  v=0 


Row  Ak(m+v)  =  0  for  ra+v+1  Sk*  ntf-v+o-1  ,  so  if  empty  suns  are  inter¬ 
preted  as  zero,  we  can  write  (2.18)  z 


c  nr+v  n+cr 

E  V'-*"''’1 '  E  *  E  W"”*  ■ 

w=K>  ^fc=0  k=mi-v+tr 


a 

a  ¥  YL  ca-w(n+v)°(®fv)  (2.19) 

v=0 


o 

a  ■  ^n+cr  *  21  co-v(B‘,‘v)n(m+'#)  a  0  . 
v*0 


For  a  »  2  ,  this  result  io  given  by  Shlntonl  (190b). 

How  suppose  that  u  of  the  R^'a  behave  similarly  as  »•*•  t  but  that 
the  ratio  of  any  one  of  these  to  each  of  the  <j-u  oilier  R,  'a  approaches 
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I 

i 

j 

i 

I 

zero  as  m-*®  ,  A  Generalization  of  a  method  due  t.o  Luke  (unpublisheu),  vhich 
was  in  turn  suceectcd  by  Clenslxav  (1964)  for  a  three-term  recursion  relation, 
can  often  be  used  to  obtain  any  one  of  the  first  u  solutions  corresponding 
to  the  u  V®  •  Clenshav  (1962)  originally  used  this  method,  for  0=2  ,  to 
compute  coefficients  for  the  Chcbyshev  polynomial  expansions  of  certain  mathe¬ 
matical  functions. 

Without  loss  of  generality,  ve  may  assume  y^(n)  is  the  solution  of 
(2.1)  that  vc  vish  to  compute.  The  algorithm  is  then  described  by  the  follow¬ 
ing  theorem. 

Theorem  2.5 

Let  the  constants  j  be  given  for  OSk<«»,l<j^u,  and  define 

m 

3  51  (2^S9,lSJiu,2«uS5. 

k=0 

(2.20) 

Let  Rjj  be  bounded  and  bounded  away  from  zero  as  m-*«»  for  2  Sh<u  , 
while  Rjf*®  °*  for  Sh  so  .  Let  rn(m)  be  as  in  (2.6)  with 

Ljj  *  Lfc^  and  lot  0(m)  ,  Ti(m)  be  non-zero  for  m  sufficiently  large.  let 
also 

lim  Sj,,JAn  *  0  * 

llM  <  *  »  2  *  h  *  u  »  *  1  '  1  *  J  *  u  *  |Vj|i  **  °  * 

1  (2.21) 
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Then  for  m^  sufficiently  largo,  vc  can  determine  ntjj  for  2  Sh  S  u  , 


<  mg  C  mj  <  . . .  <  ,  so  that  the  system  of  equations 


u  “V 

y~  Tty  ^  ^  J  s  2,3,...,U 

v=*l  k*0 


u 


H  *v  3  1  > 

v*l 


(2.23) 


has  a  unique  solution  |  rr^j  (depending,  of  course,  on  the  ) 
Let  |  Rtl(mj )  |  ^  be  bounded  away  from  aero  as  .  Then 


u 

lim  ^  nvrn(my)  =  yx(n)  ,  n  >  0  .  (2.24) 

"  v5*! 


Proof: 

Kota  that  since  ?x(my),  O(n'v)  arc  not  zero  for  a^  sufficiently  large, 
the  system  (2.23)  is  well-defined. 

Equation  (2.23)  is 


Z  *»*[  Z  M'vJSh.oKjl  -  Msis«. 

v«l  L  h«l  J 

nv*  -  ttvTi(»v^/{l(»v)  » 

The  determinant  of  (2.25)  is 

IS 


(2.25) 


(2.26) 


u  u 

(l  *  °Ui)  .  1-*  • 

The  can  be  chosen  so  that  (2.26)  is  not  zero,  else  the  are 
linearly  dependent,  by  A-I,  A- II.  Since  j  /  0  ,  ve  conclude  that 
(2.23)  uniquely  determines  for  sufficiently  large.  Furthermore, 
vlien  (2.25)  is  solved  for  ttv*  ,  one  finds  that  this  quantity  (and  hence  ny 
itself)  is  bounded  as  . 

We  write 

ti  <J 

y  TtTr^ii^*V^  ~  y  Csys(ii)  .*  cs  =  cs(®1j®2.»  *  *  •  *  (2.27) 

v=l  s=l 


Then 


Trv®s(Bv) 


-E  — 

v=l  a 
h==l 


Clearly 


(2.28) 


lim  cs=0,u+l^s^a.  (2.29) 


Also,  from  (2.25)  and  the  boundedness  of  the  rr^  we  have 

u 

y~  \  jCjj  =  1  +  o(l)  ,  Jnj-Xn  ,  IS  j  <u,  (2.30) 

h=l 
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and  since  A^  j  =  1  ,  ve  conclude  from  this  system  of  equations  that 


lim  cQ 


8  =  1  , 

1  0,  2  S  s  £  u  , 


(2.31) 


which,  when  used  with  (2.29)  in  (2.27),  gives  the  theorem. 


The  application  of  the  above  theorem  requires  that  we  know  u  normaliza¬ 
tion  relations  for  the  desired  solution,  y-^(n)  .  If,  instead,  we  know  u 
values  of  yj.  ,  then  the  following  result  can  be  used. 


Theorem  2.4 


lim  %  =  0  ;  u  +  1  £  h  *  o  , 


(2.32) 


while  is  bounded  and  bounded  away  from  zero  as  m->«  for  1  £  h  £  u  . 

Then  we  can  determine  ,  1  sh  «  u  ,  0  £  k]_  <  kg  <  k3  <  . . .  <  ku  , 
so  that 


(2.33) 


and  for  mj.  sufficiently  large,  can  be  determined,  mi  <  mg  <  ...  <  mu  , 
so  that  the  system  of  equations 


Y  =  ^(^p  »  1  *  *  u  *  (2.34) 

V=1  J  /  1 

has  a  unique  solution,  ,  1  £  h  £  u  .  Furthermore,  suppose  |  ^h(^j)  is 
bounded  away  from  zero  as  m^-*«  .  Then 


lim  Y  Mn^jA^W  =  yi(n)  >  n  *  ° 

mp-*®  v=l  / 


(2.35) 


Proof : 


As  for  Theorem  2.3. 


Unless  additional  assumptions  are  made  about  the  nature  of  the  coeffi¬ 
cients  Cv(n)  in  the  difference  equation,  it  may  not  be  possible  to  find  a 
fuidamental  set  so  that  T]_(m)  is  non-zero  for  m  large.  However,  in  most 
applied  problems,  in  particular  in  all  those  examples  of  computation  by  back¬ 
ward  recursion  which  we  discussed  in  the  introduction  of  this  work,  the  dif¬ 
ference  equation  in  question  possessed  coefficients  which  were  rational  func¬ 
tions  of  n  .  If  this  is  the  case,  as  we  shall  see,  a  fundamental  set  can 
always  be  chosen  so  that  T^(in)  is  non-zero  for  m  sufficiently  large,  and 
the  Miller  algorithm  at  least  has  a  chance  of  converging.  Even  if  it  is  only 
required  that  Cv(n)  possess  an  asymptotic  series  in  powers  of  n'1^  ,  u>  an 


IB 


,  ->-  |1 — -  -'finwi  riiii in t  " 


integer  >■  1  ,  the  r-arac  is  true. 


It  is  thus  natural  that  we  turn  our  attention  to  this  kind  of  equation. 
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III. 


In  this  and  the  following  chapters  of  this  work,  the  standard  form  for 
the  difference  equation  will  be 


a 

y~  Cv(n)y(n+v)  =  0  ,  Cc  =  1  ,  Ca  /  0  ;  o  i  2  ,  n  s  0  ,  (3.1) 

v=0 


Cv=  Cv(n)  ~  n 


{co,v 


-l/u)  -2/tu 

*  cl,vn  +  c2,vn  + 


(3.2) 


where  X^  is  an  integer,  to  is  an  integer  a  1  ,  and  cQ^v  ^  0  unless 
Cv  =  0  .  We  also  assume  that  the  coefficients  sure  written  with  the  smallest 
possible  value  of  u>  .  Note  that  the  equation  adjoint  to  (3.1)  also  has  co¬ 
efficients  of  the  same  general  form  as  (3.2).  2y  (A. 12),  equation  (3.1)  has 
a  unique  solution  y(n)  satisfying 


y(nQ+j)  =ffj,lsj^a,n0^0  .  (3.3) 

Thus,  for  n  a  nQ  2  0  ,  (3.1)  has  a  fundamental  set. 

We  shall  see  that  there  exists  a  certain  fundamental  set  for  (3.1) 
whose  members  share  an  unusual  property:  each  has  an  asymptotic  expansion, 
valid  as  n— ,  which  consists  of  an  exponential  leading  term  multiplied  by 
a  descending  series  of  the  kind  (3.2)  (where,  however,  m  may  be  replaced  by 
an  integral  multiple  of  <u  ,  see  (3.4)-(3.7)  below).  Essentially,  these 
series  are  the  same  as  the  so-called  subnormal  series  encountered  in  the  study 
of  ordinary  linear  differential  equations  with  polynomial  coefficients  near 
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rJ>;  ‘ : 


singular  points,  see  Ince  (1956,  Ch.  1?).  This  is  another  example  of  the 
close  analogy  between  differential  equations  and  difference  equations. 

For  our  purposes,  the  existence  of  solutions  of  the  difference  equation 
(3.1)  which  have  such  asymptotic  representations  is  quite  important,  since  the 
very  form  of  the  asymptotic  series  enables  us  to  determine  much  more  practical 
conditions  for  the  convergence  of  the  class  of  algorithms  discussed  in 
Chapter  II. 

However,  we  will  have  to  examine  the  properties,  algebraic  and  analytic, 
of  these  solutions  in  great  detail  before  we  can  attack  the  convergence  prob¬ 
lem  directly,  and  the  present  chapter  is  devoted  entirely  to  this  study. 

YJe  begin  with  several  definitions. 

Consider  the  series 


r 

Q(n)  =  Q(pjn)  =  u0n  In  n  +  p-jn 


P+1-3 


(3.4) 

(3.5) 


s(n)=  s(p;n)  =  n9  ^  (in  n )Jn  q.(p;n)  , 

j=0 


(3.6) 


<lj(p;n)  =  qj(n)  =  ">_  bSjJn  S//p  , 
s=0 


(3.7) 


where  p  ,  r^  ,  pQp  are  integers,  p  £  1  ,  Uj  ,  0  ,  bg  j  are  complex  param- 
eters,  b0^j  /  0  ,  unless  bg^j  ~  0  for  all  0  £  s  <  »  ,  ro®0, 


-tt  £  Im  ui  <  tt  • 
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Definition  3.1 

The  series  (3.4),  called  a  formal  series  (F.S.),  will  he  called  a  formal 
scries  solution  (F.S.S. )  of  (3.1)  if,  when  it  is  substituted  in  (3.1),  the 

Q  ^  ri  \  . 

equation  is  divided  by  e  '  7  and  the  obvious  algebraic  manipulations  (see 
below)  are  performed,  then  the  coefficient  of  each  quantity 


n<*rW«(ln  n)j 


0  s  j  s  t  ,  r,s  -  0,tl,±2, 


(3.8) 


is  equal  to  zero. 


A  concept  of  formal  equality  between  two  F.S.  can  be  defined  by  requiring 
that,  when  the  series  are  written  with  the  same  value  of  p  (as  is  always  pos¬ 
sible),  then  the  parameters  t  ,  0  ,  bs^j  ,  r^  ,  for  both  series  are  the 
same.  Formal  equality  of  formal  solutions  also  ar'ses  in  the  theory  of  or¬ 
dinary  differential  equations  (see  Coddington  and  Levinson  (1955),  p.  114  ff). 

The  construction  of  F.S.S.  may  be  carried  out  by  using  the  identities 


u(n+v)  _ 


X 

h,J=l 


ak,jn 


k(6-l)+l-j 


(3.9) 


(3.10) 

(3.11) 
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[ln(n+v)J  = 


In  n  +  — 
n 


2rf 


(3.12) 


although,  in  practice,  it  is  difficult  to  obtain  by  hand  other  than  the  first 
few  terms  this  way,  see  Birkhoff  (1930).* 

Very  often,  wc  shall  let  "Q(n),"  "s(n),"  be  generic  symbols  for  the 
series  (3.5),  (3.6).  The  series  so  denoted  do  not  necessarily  have  the  same 
values  of  the  parameters  9  ,  t  ,  bs  ^  ,  pj  ,  when  the  series  occur  in  dif¬ 
ferent  equations.  If,  however,  it  is  necessary  to  differentiate  between  two 
such  series,  we  shall  do  so  with  subscripts,  e.g.,  Q^(n)  ,  s^(n)  ,  Qg(n)  , 
s2(n),  etc.  With  this  convention  understood,  we  see  that  F.S.  possess  the 
following  properties: 

e°'(n+V)s^(n+v)  =  e^n^s2(n)  ,  v  =  1,2,...,  (3.13) 


QL(pjn)  Q2(p';n)  Q3(p*;n) 

e  s-L(p;n)  •  e  s2(p';n)  =  e  s3(p*;n)  ,  (3.14) 


where  p*  is  the  least  common  multiple  of  p  and  p '  . 

The  sum  of  two  F.S.  is  not  in  general  a  F.S.,  but  if  Q(n)  ,  0  are  the 
same  .for  both  series,  we  have 


*  J.C.P.  Miller  has  brought  to  our  attention  the  fact  that  the  determination 
of  these  series  can  be  done  very  efficiently  by  computers. 


23 


(3.15) 


eQ(n)n9q1(p;n)  +  eQ^n^n0qg(p' jn)  =>  e^n^n9q3(p*;n)  . 


Definition  3.2 


f(n)  ~  e^n^s(n)  ,  n—»«»  , 


(3.16) 


means  that,  for  every  k  a  1  ,  we  can  determine  functions  j(n)  , 
0  s  j  st,  such  that 

e^K-V)  =  f_  (In  n)Jnrt-J/P  £  b, 

po  s=o  * 


+  n  ^  (in  n)^n  ^  ^  Ak  .(n) 

j~0 


(3.17) 


and  lAk  j|  is  bounded  as  n— ,  for  all  k,j  . 


See  Birkhoff  and  Trjitzinsky  (1932,  p.  62).  If  t  -  0  ,  this  definition 
coincides  with  (A.7)-(A.8).  Also  (3.17)  is  unique,  since  it  is  readily  veri¬ 
fied  that  zero  has  no  non-trivial  representation  of  the  form  (3.16). 

Definition  3.3 


Let 


Wk  = 


Qh(n+J  )  /  , . v 

e  sh(n+j) 


(3.18) 
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By  (3.12)-(3.14),  Wk 


is  a  F.S.,  and 


Wk  *>  exp 


XV 

Z 

j=i 


9j(») 


\  s(n)  «  e^nV{a) 


(3.19) 


•  S.  | 


Qh(n) 


We  say  the  k  F.S.  se  s.  (n) 


are  formally  linearly  independ<  nt  if 


s(n)  ^  0  .  Otherwise,  they  are  formally  linearly  dependent. 


Definition  3.4 

There  exist  exactly  r  F.S. 3.  of  a  certain  type  (e.g.,  with  Q(n)  s  0) 
for  the  equation  (3.1)  means  r  F.S.  of  that  type  can  be  constructed  which  are 
formally  linearly  independent,  and  any  r+1  such  F.S.S.  are  formally  linearly 
dependent. 


Now  we  may  formulate  two  very  important  questions  about  the  difference 
equation  (3.1).  Does  the  equation  always  possess  exactly  o  F.S.S.  of  the 
general  type  (3.4)?  If  so,  what  relationship  do  these  F.S.S.  bear  to  the 
fundamental  sets  for  the  equation? 

These  questions  were  answered  partially  by  a  number  of  mathematicians, 
see  Adams  (1928)  and  the  references  given  there.  Hcvever,  only  with  the  ad¬ 
vent  of  two  papers,  the  first  by  Birkhoff  (1930)  and  the  second  by  Birkhoff 
and  Trjitzlnsky  (1932)  was  the  theory  completed.  The  results  of  these  two 

2? 


writers  yield 


Theorem  3.1  (Blrkhoff-TrJitKinsky) 

There  exist  exactly  o  F.S.S.  of  equation  (3.1)  of  type  (3.4),  where 
p  ■  **u>  ,  for  some  integer  V*  1  ,  and  each  F.S.S.  represents  asymptotically 
some  solution  of  the  equation  in  the  sense  of  (3.16).  Ifce  o  solutions  so 
represented  constitute  a  fundamental  set  for  the  equation. 


Definition  3.5 

The  particular  fundamental  sets  mentioned  in  Theorem  3.1  will  be  called 
Blrkhoff  sets.  Each  member  of  a  Birkhoff  set  is  called  a  Birkhoff  solution. 


Now  let  {yh(n)J  ,  1  s  h  s  o  ,  be  a  Birkhoff  set  for  (3.1),  i.e.,  let 


Qh(p;n) 

yh(n)  ~  e  sh(p;r  )  ,  n— >*  ,  1  s  h  s  a  , 


(3.20) 


which  is  permissible,  since  we  can  write  all  the  F.S.S.  of  Theorem  3.1  with  a 
common  value  of  p  . 


Then  none  of  the  determinants 


-  yBh(n+<} )  ,  1  «  <  Wg  <  . . .  <  *  o  , 


(3.21) 


can  be  zero  for  n  sufficiently  large,  and  in  particular,  y^(n)  i  0  far  n 
sufficiently  large,  1  i  h  i  a  .  In  fact,  since  the  F.S.S.  are  formally 


linearly  independent, 


is  of  the  fora  m/p  ,  m  an  integer,  so  the  difference  equation  for 
y(n)  will  be  c*  the  same  kind  as  (3.1),  except  with  w  replaced  by  p  .  Let 
Cyjn)  be  the  coefficient  of  the  transformed  equation  and 


cp^Cn)  *  n9  21  (In  n)iqt.^_k(p;n)nrk“^P(t+l-k )jJ  l'.  ,  (3.86) 

0 

<pt(n)  being,  by  hypothesis,  a  F.S.S.  of  that  equation.  Substituting  (3.26) 
with  k  -  t  into  (3.1)  and  setting  to  zero  the  coefficient  of  (in  n)U  gives 


2~  Cv(n)(n+v)9  21  qjt+u ( P J n+v ) (u+ 1  )x ( n+v )  t  X'UL(l+v/n)l  /  ll  =  0, 

v=0  4=0 


(3.27) 


0  £  u  £  t 


Let  u  =  t-j  ; 

22  C^nHn+v)9  2Z  q4+t-j(n+v)(n+v)  J  1  (t+l-j)x[ln(l+v/n)j  / Ll  =0 
v=0  4=0  J 

OSj  St  . 

Nov  writing 


(3.28) 


£in(l>-v/n)j  =  J\ui(n+v)  -  In  nj 


(3.29) 


and  expanding  by  the  binomial  theorem,  we  get 
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_  a  i  „  J-u  p"  r,  ... 

21  Cv(n)(n+v)  V  (In  n)  (-)  /u.'  Y"  <lt+4+u.j(n+v)(n+v) 
v=0  u=£)  i=0 


X  (t.+1-j  )A+u[m(n+v)]A/j«  *o,osjst  . 


From  (3.26)  we  have 


(3.30) 


J-u 


lr 

<Pj_u(n)  -  n9  21  (In  Q)\+i+u.j(n)n  ^^(t+l+u-jj/li  ,  (3.31) 

t~  0 


and  so 


a  J 

21  cv(n)  27  ^'^ln  n)U(t+l^)u9j_u(n+v)/u!  =0  ,  Osjst.  (3.32) 
v=0  u=-0 


Now  suppose  cp^n^cp^n), . .  .,9^(0)  are  F.S.S.  of  the  transformed  equa¬ 
tion,  and  let  j  =  k+1  in  (3.32).  We  have  \ 


X.  Cv(n)^+1(n+v)  =  0 
v=0 


(3.33) 


But  for  j  -  0  ,  (3.28)  gives 


0  .  a  _ 

21  Cv(n)(n+v)  ^(n+v)  =  27  cv(nWn+v^  =  0  » 

v=0  v=0 


(3.34) 


so  the  induction  is  complete,  and  the  ^(n)  ,  0  s  h  s  t  ,  satisfy  (3.1). 


“*  ^«**nw«f**aV*~vssrii  2 
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Birkhoff  (  J30)  has  noted  that,  once  a  maximum  value  of  t  is  found  so 

(L^n) 

that  e^  sQ(n)  is  a  F.S.S.  of  (3.1),  then  there  are  t  other  F.S.S., 
Q_(n) 

e  8^(0)  >  1  £  h  s  t  ,  and  in  sh(n)  ,  In  n  occurs  to  the  t-h  power. 
However,  the  explicit  form  of  the  series,  i.e.,  (3.23),  seems  to  be  new. 

We  will  need  information  about  the  number  and  form  of  those  F.S.S.  of 
(3.1)  whose  exponential  leading  terms  are  constant.  The  necessary  results 
are  contained  in 

Theorem  3.3 

Let  (3.l)-(3.2)  hold,  and  define 


Pk~  pk(n)  *  X!  >  k  =  1>2'***' 

v=l 


po  3  »<.(■>)- II  e.  • 

v=0 


(3.55) 


Then  ve  can  write 


Pk  ~  nafc/U)[a0>k  +  al#R  n1^  +  +  •••],  n-»»  ,  (3.36) 


»k  an  integer,  aQ  k  f  0  unless  Pk  =  0  ,  in  which  case  ve  interpret 
ak  =  *  ***• 
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t  =»  max  (o'k/o>  -  k) 
0  k  k  so 


(3.37) 


and  let 

kQ  <  k^  <  kQ  <  •  •  •  <  kj  (3.38) 

be  those  values  of  k  for  which 

-  k  =  T  ,  Osk<«  .  (3.39) 

Then  s  a  ,  and  there  exist  exactly  k^  F.S.S.  of  (3.1)  of  the  form 
s(p;n)  ,  and  p  =  «  .  Each  represents  asymptotically  a  solution  of  (3.1)  as 

n — ,  the  k.  solutions  so  represented  being  linearly  independent. 

? 

Furthermore,  0  =  8h  is  one  of  the  k^  values  satisfying 

♦  k. 

Go(0)  -  (*)  H^^k/V  "  0  *  (3.40) 

and  if  no  two  of  the  0h's  differ  by  integral  multiples  of  l/uo  ,  then  the 
F.S.S.  take  the  form 

s(n)  =  n0[bo  +  bjn” L^<B  +...],  b0  t  0  .  (3.41) 


Logarithmic  solutions  can  occur  only  if  seme  of  the  roots  of  Go(0)  differ 
by  integral  multiples  of  l/u>  . 


Preview  of  Proof: 


First,  ve  shew  that  p  =  <u  in  any  purely  algebraic -logarithmic  F.S.S. 
of  the  difference  equation  (3.1). 

Second,  by  considering  (3.35)  as  a  system  of  equations  in  the  "unknowns" 
Cy(n)  ,  0  s  »  s  o  ,  ve  shew  there  exists  a  finite  t  satisfying  (3.37)  and 
that  kj  £  a  . 

Third,  the  existence  of  F.S.S.  is  shown  by  actual  construction,  using 
Frobenius’  method,  first  for  the  case  where  none  of  the  roots  of  Go(0 ) 
differ  by  integral  multiples  of  l/u)  ,  and  next  far  the  case  where  none  of 
the  roots  of  Go(0)  are  equal,  but  where  there  exists  a  subset  of  roots 
differing  from  each  other  by  integral  multiplies  of  l/a>  .  It  is  then  shown 
that  in  these  two  cases,  any  purely  algebraic -logarithmic  F.S.S.  of  (3.1)  can 
be  expressed  as  a  linear  combination  of  the  solutions  already  constructed, 
i.e.,  for  these  two  cases,  there  are  exactly  k^  such  F.S.S. 

lastly,  we  indicate  briefly  how  similar  results  are  obtained  when  some 
of  the  roots  of  GQ(9)  are  equal,  and  the  solutions  for  this  case  are 
displayed. 


Proof: 

Let 


A8  jr  (1A  A^AhPiA) 

x=o 


(3.42) 
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The  difference  equation  is  unchanged  by  replacing  n  by  ne2^nitu  ,  h 
an  integer,  so  each  of  the  functions  cp^[ne^^n^UJ]  is  also  a  F.S.S.  of  (3.1). 
An  application  of  Theorem  3.2,  or  a  result  of  Birkhoff  (1930,  section  2), 
shows,  furthermore,  that  each  of  the  functions 


Vk(a)  =  X  e2h"1J /*'  sjA(“in) 

j=0 


(3.48) 


is  a  F.S.S.  of  (3.1).  V/e  can  determine  unique  constants  A^j  ,  O^h^v-1  , 
where  v  is  the  number  of  non-zero  Sj  ^  ,  so  that 

v-1 

51  Ah,A,k(n)  3  sj,k(tu;n)  >  0  s  j  s  v-1  (3.49) 

h^O 


since  the  determinant  of  the  system  of  equations  (3.48)  is  a  Vandermonde 
determinant  which  is  non- zero. 

But  (3.49),  when  substituted  in  (3.45),  shows  that  every  purely  algebraic 
logarithmic  F.S.S.  of  (3.1)  can  be  expressed  as  a  finite  linear  combintation 
of  similar  solutions,  each  of'  which  can  be  written  with  p  -  uu  ,  and  the 
first  step  of  the  proof  is  completed. 

Not  all  the  P^'s  ,  0  s  k  s  o  ,  can  be  zero  unless  CQ  =  C-^  -  ...  = 

C  »  0  ,  since  the  determinant 
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does  not  vanish.  Hence  there  exists  at  least  one  finite  ,  0  £  k  s  o  , 
and  therefore  a  finite  r  satisfying  (3.37).  Also,  from  (3.35),,  we  have 
for  k  =  1,2, , 


and 


m; 


P  12 
o+k  L 


a+k 


a+k 


=  0 


(3.51) 


Pa+k  "  *l,kPl  +  i2,kP2+  •**  +  1q, kpa  ' 


(3.52) 


ar^/tu  s  max  cr./w  ,  k  =  1,2,..., 
1  s  s  a  J 


(3.53) 


or 


<VkAu  ~  (o+k) 

-  * 


<  max 
1  «j  Sd 


[«j/w  -  J  >  k  -  1,2,..., 


(3.54) 


and  k.  s  a  . 

♦ 

Now  assume 


y(n) 


9  V- 


n 


Zv 

s=0 


■s/i 


to 


bo  t  0 


(3.55) 
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(3.56) 


y(n+v)  ~n°  bsn'S//w  Y  (s/oM))kHkvkA!nk 


s=0  k=0 


Substituting  (3.56)  in  (3.1)  gives 


03  ,  CO 

0  ST“  ,  -s/< 


n9  Y  V’S/U)  21  (sA-8)k(-)kPk/k.‘nk  =  0  .  (3.57) 


s=0  k=0 


We  can  write 


«  00 

21  (s/u)~9)k(-)  PkA-'nk  =  nT  21  Gje-s/uOn"1"/10  ,  (3.58) 


k=0 


m=0 


and  so 


a  .  m 

2]  Cvy(n+v)  ~  n  T  2.  n_m Y  bsGm-s^~s^  =  0  *  (3.59) 


v=0 


m=0  s=0 


We  must  have 


m 

s 


^s^m-s (9~® A )  —  0  ,  m  —  0, 1*2, . . .  . 


(3.60) 


If  none  of  the  roots,  0h  ,  of  Go(0)  differ  by  integral  multiples  of  l/uj  , 
then  the  construction  of  the  k,^  solutions  can  proceed  directly  from  (3. 60). 
If  some  of  the  roots  of  GQ(9 )  do  differ  by  integral  multiples  of  l/uj  ,  then 
the  construction  of  solutions  is  done  by  the  method  of  Frobenius,  as  follows. 
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We  now  differentiate  (3.69)  h  times  with  respect  to  9  ,  set  0  =  0h+lJ 


and  use  the  fact  that 


^1.(9  )l 


=  0  ,  0  s  s  s:  > 


(3.72) 


to  obtain  the  X+l  F.S.S. 


9.  JL 


<Ph+L(n)  =  n  1  X  (-h)x(-/(ln  n)h  lgJ,  h(n)n  V  / ,  (3.73) 

1=0 


CO  ^ 

=  2_  bs,.i,hn  S  U  *  bs,i,h  =  “ ~l 

s=0  S9*  ^  1  9=9. 


0  s  h  s:  x 


(3.74) 


Since 


'.h  _ 

Vh,h  =  @o(0) 

J  '  agh 


t  o  , 


(3.75) 


the  above  solutions  are  clearly  formally  linearly  independent,  by  A-VII. 

Let  (3.42),  with  p  =  uj  ,  be  a  F.S.S.  of  (3.1  and  let  GQ(9)  contain 
no  multiple  roots.  We  wish  to  show  tp^n)  can  be  expressed  as  a  linear  com¬ 
bination  of  the  solutions  already  constructed.  Denote  by  ds^  the  coeffi¬ 
cient  of  n  s ^  in  q^(n)  .  We  can  write 


I 


(3.76) 


$k(n)  =  21  ^7  (»3) X  t/ 

C=0  99x 

and  substitute  this  in  (3.1),  we  find  that  9  and  4g  j  must  be  such  that 


0  £  v  £  k  •  (3.7?) 

For  v  =  0  ,  the  above  equation  demands  that  Gc(9 )  -  0  .  Thus  assume 
9  =  0X  ,  since  the  coses  6  =  9h  2  s  h  £  x+1  ,  can  be  treated  similarly, 
as  can  th^  cases  where  8  belongs  to  any  other  group  of  roots  of  G0(e) 
which  differ  by  integral  multiples  of  l/m  .  If  8  =  9  .  where  0  differs 

|X  i* 

from  no  uther  root  of  GQ(3)  by  an  integral  multiple  of  l/oo  ,  then  the 
analysis  is  quite  simple;  we  must  have  k  =  0  in  (3.76),  and  the  result  is 
one  of  the  purely  algebraic  F.S.S.  (3 .41). 

One  can  show  by  induction,  using  formula  (3.77),  that  the  r-a's  cannot 
ue  arbitrary,  but  must  satisfy  the  relations 

0  =  rn  =  L,  <  r-.  -  Lr  <  ro  =  Lj  <  ...  <  rk  =  L,  s  L .  ,  (3.78) 

and  <)£<•••  .  (To  avoid  double  subscripts  on  the  r^'s  ,  it 

is  assumed  that  any  vhich  corresponds  to  a  <iyk-j1(n)  =  0  is  deleted  from 
the  chain  (3.78).  ) 
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How  let  |»  be  the  highest  power  or  In  n  which  occurs  in  any  of  the 


functions  ^(n)  ,  1  *  h  s  x+1  -  Hote  t  iat  the  last  term  in  is 

(3.79) 


0/k,/® 

n  «h,h<n>  '  8h,h<°>  t  0  * 


Hence  we  may  select  a  (^(n)  which  contains  the  term  n  A  K  gy  y(n)  , 

and  subtracting  a  constant  multiple  of  this  function  from  9fc(n)  leaves  a 

8.+r*/ta  * 

series  of  the  sane  kind  whose  lead  torn  is  either  zero  or  n  L  *  q0(n)  , 

rk  <  rk  *  But>  by  (3.78),  r£  oust  be  one  of  the  1^  .  Hence  we  can  find 

8,+r* 

another  cp^n)  with  lead  term  n  i  *  and  subtract  a  constant  multiple  of 

this  series  from  the  above  series  so  that  the  initial  power  of  n  in  the 

9ihJ*  _  -  . 

lead  term  is  again  diminished,  to  n  >  rk  <  rk  *  be  con¬ 

tinued  until  ^(n)  is  reduced  to 


1 


^(n)  =  fk(n)  +  gk(n) 


(3.80) 


where  gk(n)  is  a  linear  combination  of  the  ^(n)  and 


fk(n)  =  n  L  (In  n)\j(«M;n)n 
i=l 


x  ?k -V® 

•n  !r»  * 


k  =  max 


(3.81) 


also  satisfies  the  difference  equation.  It  is  our  intention  to  prove  that 
any  F.S.S.  of  the  form  (3.81)  must  be  identically  zero. 

Denote  by  d_  ,  the  coefficient  of  n  S^U)  in  q.(n)  .  Since  d„  =  0  , 

S y  A?  X  O 

the  result  of  writing  out  (3.77)  for  (3.81)  and  k  =  v  =  k  is 
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A 

. — 


-U »■*»*$*  tfCrivVfjp*’'* 


I  "V“  t  ”W”  t  Vt,  ^  Vs  (»  * 

X=0  m=0  5=0  V  *"  / 


0=0, 


=  0  ,  (3.82) 


rr  ,/u) 

and  setting  to  zero  the  coefficient  of  the  highest  power  of  n  (n  K  ) 
gives 


Vi  J;  =  0 


(3.83) 


which  is  only  possible,  since  GQ(0 )  has  no  double  roots,  if  dQ  ^  =  0  . 

Ifence,  by  (3-7),  q^(n)  =  0  .  The  suai  (3.81)  becomes  a  sum  from  l  =  2  to  k. 
New  write  out  (3.77)  for  v  =  k-1  =  k-1  to  show  that  qg(n)  S  0  .  Eventually, 
we  arrive  at  q^n)  =  qg(n)  =  ...  =  q^(n)  =  0  ,  or  fk(n)  =  0  ,  so  ^(n) 
is  a  linear  combination  cf  the  cp^n)  ,  1  i  h  s  X+i  >  and  hence  k  s  p,  . 

Thus,  if  none  of  the  roots  cf  Gft(6)  are  eoual,  there  exist  exactly  k 

7 

algebraic  logarithmic  solutions  of  (3.1). 

The  same  analysis  can  be  conducted  when  some  of  the  9-n's  are  allowed 
to  be  equal,  and  the  construction  of  the  solutions  in  this  case  again  is 
analogous  to  the  procedure  used  for  differential  equations,  see  Ince  (1956, 

Ch-  XVI)  or  Forsyth  (1902,  Ch.  II).  Here  we  simply  display  the  solutions, 
since  their  forms  will  be  of  importance  later. 

Let  0^  be  as  in  (3.65)-(3.66 ),  and,  furthermore,  let  it  be  a  root  of 
Gq(9)  of  multiplicity  6^  .'i  1  .  Define 


42 


(3.84) 


6h  =  6L  +  62  +  •••  +  6h  ,  6*  =  °  . 


The  6jj  F.S.S.  corresponding  to  0^  axe 


CPh,j_6h-l+1 


/  \  9h  'T-  ("J  ,,  vj-v  V-  -s/iu  dV  —  t*\ 

(n)  =  n  n  >  - -  (In  n)  Z_n  7  3S(0)  ' 

v=0  *  s=0  90  9=9h 


5h-i  s  J  s  * 


(3.85) 


Note  that 


J  - 


3O(0) 


(3.86) 


since  the  multiplicity  of  (9-0^)  in  gQ(0)  is  exactly  6^^  - 

The  linear  independence  of  the  above  solutions  can  be  shown  by  the  same 
arguments  used  by  Ince,  p.  402,  and  the  proof  that  there  are  always  exactly 
kj  such  F.S.S.  requires  the  same  type  of  reasoning  as  when  the  roots  of 
GQ(0)  are  simple,  but  the  details  are  considerably  more  messy. 

Theorem  3.1  is  now  applied  to  show  that  the  F.S.S.  constructed  above 
represent  asymptotically  k^  linearly  independent  solutions  of  (3.1),  and 
the  proof  of  Theorem  3.3  is  complete. 


Definition  3.6 


A  Birkhoff  set  {yh(n)j  for  equation  (3.1)  is  a  canonical  set  if 
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t .  t 


where  y^n)  is  any  member  of  3  other  than  z^  and  y2(n)  is  any  member 
of  S  other  than  z2  . 

For  the  construction  (3.85),  z^  and  z2  correspond  to  F.S.S.  ^^(n) 

and  9^+i  g  (n)  ,  respectively,  and 
X+1 


(1)]  , 

(3.91) 

[i  +  o(i)]  . 

(3.92) 

z2(n)  =  c2n  X+l(ln  n)  X+1 

(in  the  general  case,  a  term  e^n^  will  appear  on  the  right-hand  side 
of  (3.9l)-(3.92). ) 

Theorem  3.4 


Let  {yh(n)|  be  a  Birkhoff  set  for  (3.1), 


yh(n)  ~  e  n  sh(p;n)  ,  1  £  h  s  o 


(3.93) 


Then: 

i) 

where  in  Q(n) 


D(n)  ~  e^^u),n^s(w;n)  , 


(3.94) 


“o  =  -K>  '  ^  •"  <->°e  ^  /co,o  '  ^2  ' 


U) 

^3  -7~ 


U-2) 


[(cl,a/2co,o)  "  (c2,a/co,a)] 


(3.95) 
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where  Q  is  as  in  (3.2),  and  s(to;n)  is  free  from  logarithms; 


CT 


ii) 

Q(u>;n)  =  Yi  ^h^p;n^  * 

h=l 

(3\  96  ) 

iii) 

_  -Qh(n) 

yh(n)  =  Th(n)/D(n)  ~  e  sh(p;n)  ,  lshsff, 

(3.97) 

are  a  set  of  Birkhoff  solutions  for  the  equation  adjoint  to  (3.1). 


Proof: 

The  proofs  of  i)  and  iii)  are  purely  computational,  i)  following  from 
the  difference  equation  for  D(n)  ,  (A. 15).  No  logarithms  appear  in  s(u>;n) 
since,  by  Theorem  3.2,  the  F.S.S.  of  a  first  order  difference  equation  can 
never  contain  logarithms. 

To  prove  iii),  note  that  from  A-VI, 

frh(n+j - 1  )/D(n+j -1  )J ^ X  Jyj (n+h-1  )1  ^  =  fe jh| ^  >  (3.98) 


where 


c 


Jh 


fl 


>  J  -  h  , 


> 


0  ,  h  >  j 


(3.99) 


and  so 


(3.100) 


Since  the  y^(n)  can  *>e  represented  by  the  F.S.  (3.97)  which  are 
formally  linearly  independent,  by  (3.100),  it  follows  that  they  are  Birkhoff 
solutions. 


We  close  this  chapter  with  a  theorem  on  exponential  sums. 


Theorem  3.5 
Let 

n 

S  =  S(n)  =  21  f(k)  *  n  2  0  >  (3.101) 

k=0 


where 


f(n)  =  e^n^v(n)[l  +  o(l)]  ,  n— >«  ,  (3.102) 

v(n)  =  n®(ln  n)r  ,  (3.106) 


r  is  a  non-negative  integer,  and  Q(n)  is  given  by  (3.5). 

Let  p*  be  the  first  non-zero  element  in  the  sequence  |Re  p^j  > 
h  =  0, 1,2, . . . ,p  and  p*  =  0  if  all  the  Re  ph  are  zero. 
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Then : 


i)  if  n*  <  0  ,  we  can  write 


A  -  Z  *oo  > 

k=0 


S(n)  «  A  +  0  jeQ(n+1)nv(n) J  , 


n— ; 


ii)  if  p*  >  0  , 


S(n)  =  0  je^n^nv(n)  j  , 


n— >»  ; 


iii )  if  p*  =  0  , 


0  |nv(n)j  ,  Re  9  >  -1  ; 

0  [(In  n)r+1'|  ,  Re  6  =  -1  ; 
0(1)  ,  Re  0  <  -1  . 


(3.104) 


(3.105) 


(3.106) 


(3.107) 


Proof : 


We  can  write 


f(n)  =  ce^(n)v(n)  [l  +  K(n)J  ,  n  a  1  , 


(3.10Q) 


where  K(n)— - >0  as  n— >»  . 


%\  •  5s .  -¥>; 


For  n*  <  0  )  the  series  converges  absolutely.  Define 


CO 

T(n)  =  T  f  (k)  .  (3.109) 

k=n+l 

We  have 


|e'Q(n+l)T(n)|  S  |cj  g(k)k“V|l  +  K(k)|  (3.110) 

k=n+l 


sc'  g(k)k"v  ,  (3.111) 

k=n+l 


where 


K)  =  |v(k)|  kV=Re  Q(k)"Re  Q(n+1) 


(3.112) 


and  we  take  v  real  and  >  1  . 

Now 

=  g(x)  [Re  Q*(x)  +  +  0  +  v)  1  .  (3.113) 

dx  (  x  \  In  x  J  J 

But  g(x)  f  0  for  x  >  1  ,  and  any  zero  of  the  quantity  in  brackets 
above  cannot  depend  on  n  .  Since  g(x) — >0  as  x — >»  ,  it  follows  that  g 
is  monotone  decreasing,  x  >  xQ  ,  or,  for  n  sufficiently  large 

g(k)  s  (n+l)V| v(n+l)|  s  MnV|v(n)|  ,  k  £  n+1  ,  (3.114) 
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and  so 


a-®(n+1)I(n)|  s  M  | v(n) |  nv  k'v 

k=n+l 


(3.115) 


kS|v(„)|  nv  /  x^dx  -  Ssdsi  (  (3.116) 

i  (v-L) 


which  gives  i). 

The  result  ii)  follows  by  a  simple  majorization,  while  for  iii)  it 
suffices  to  consider 


n  n 

2jv(k)|  *  (In  n)r  21  k  .  (3.117) 

k=l  k=l 


When  Re  0  z  0  ,  the  series  on  the  right  is  easily  bounded.  Otherwise 
we  use 


2  k^e  ®  <  1  +  J  xPe  ®dx  ,  Re  0  <  0  .  (3.118) 

k=l  1 
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IV. 


By  using  the  asymptotic  theory  of  the  last  chapter,  we  can  now  apply  the 
general  theory  of  computation  by  backward  recursion,  as  developed  in  Chapter 
II,  to  difference  equations  of  the  kind  (3.1). 

Throughout  this  chapter,  we  assume  the  difference  equation  in  question 
to  be  (3.1). 

Theorem  4.1 

There  exists  a  canonical  set  |cp^(n)|  for  the  equation  adjoint  to  (3.1) 
and  an  integer  u  s  1  such  that 

fch  /  0  ,  1  s  h  £u  , 

lim  cp*(m)/<p*(m)|  =  <  (4.1) 

m->»  1  [o,  u+1  shsff  . 

Furthermore,  if  u  =  1  then  for  some  n*  ,  0  s  n*  £  a-1  , 

lim  An(m)/An*(m)  =  y(n)  ,  n  s  n*  ,  (4.2) 

exists,  is  not  identically  zero  and  satisfies  (3.1),  while,  if  u  =  2  ,  then 
ml  *  ^  cau  ch°sen  50  that  (2.35)  holds  for  some  solution,  y^(n)  ,  of 
(3.1).  Furthermore,  y-j_(n)  ,  the  function  to  which  the  algorithm  converges, 
is  independent  of  the  particular  set  <p£(n)  . 
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Proof; 


Equation  (4.1)  is  obvious,  since,  clearly,  the  absolute  value  of  the 
ratio  of  any  two  canonical  solutions  must  either  approach  zero,  a  constant, 
or  infinity.  Let  {y^n)}  be  that  (unique)  fundamental  set  for  the  original 
equation  (3.1)  which  is  obtained  by  iettin*.  Th(n)/D(n)  ■  9^(n)  in  (A. 18), 
then  solving  for  y^(n )  for  n  a  a-1  ,  and  finally  usiig  (3.1)  to  compute 
yh(n)  for  all  n  2  0  ,  1  £  h  so.  We  have 


o 

An(m)  *  Y  c?h(m)yh(n) 
h=I 


(4.3) 


euid  the  statement  for  u  *  1  follows  immediately.  Note  that  An*(m)  cannot 
be  ze^o  for  more  than  a-1  consecutive  values  of  n*,  for  it  satisfies  (3.1) 
and  is  not  identically  zero.  Since  An(m)  is  fully  determined  by  the  condi¬ 
tions  (2.2),  it  is  independent  of  <pj*(n)  ,  and  so  is  yj_(n)  • 


.  i2  1  crei^rai^0(l  +  o(l)) 


(4.4) 


where  Q  ,  8  are  real,  [orj  a  c^  ,  and  in  Q  is  zero.  Alro  we  may 
assume  -n  s  ^  <  rr  ,  or  else  replace  iQ  by  iQ  i  2mi  ,  r  a  suitably 
chosen  integer. 
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One  finds 


and 


|^(m1) |  =  2cl 


sin  s 


- 2 -  +  |  (ln  m2,'ln  ml) 


+  o(l)  .  (4.5) 


First,  assuaie  Q(m)  ^  const. ,  »  0  ,  so 


1-r 

p  i  +  ±-£ 

Q(m)  =  2]p-rm  P  ,2sksp, 
r=k 


(4.6) 


Let 


mg 


K 


r  6ik 

--  +  Km^p  =  +  Kmj_^  +  0(1)  , 

*  h/C^  (i  *  ¥))  |  < 


>  (4.7) 


[xj  means  the  largest  integer  not  greater  than  x  .  We  have 


QOngHK)  i+-r  1+ 0 

— i — *£*{*  -■! 

a  &  2^  (l  +  nr) “7"  C1  *  °<mlvp>) 

-  jr  K  (l  +  tO  +  0(ail/p)  -  i  |  +  0(mIL/p)  , 


(4.8) 
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£10 


lim  | A(m^)  |  »  2c^  , 


(4.9) 


mi — >  °» 


which  assures  us  that  (2.34)  has  a  unique  solution  tT]_,  and  also  that 
(2.35)  holds. 

Next,  we  note  that  if  Q  is  constant,  we  cannot  have  9=0,  else  cp* 
is  proportional  to  <p*  .  Thus,  let 


."/lei  + 


o(D  . 


(4.10) 


Hien  (4.9)  again  holds. 


Lastly,  if  nj.  t  0  in  Q  >  let  nig  =  m^+r,  r  a  positive  integer. 
|4(m].)  |  *  ci  | e^1  -  1  +  o(l)  |  . 


(4-U) 


irpi  .  .  . 

Since  r  can  be  chosen  so  that  e  ‘  f  l  ,  (4.9)  again  follows,  and. 


by  Theorem  2.4,  the  proof  is  complete. 


Theorem  4.2 


Let  there  exist  a  canonical  set  jyh(n)|  for  (3*1)  such  that  one  of  its 
members,  say  y^(n)  ,  has  the  following  property: 


ft**/  V&Q^-W*-*****  if*1 


; 


(4.12) 


lim  y^(n)nX/yh(n)  a  0,  2  shsij 
n  — ♦» 


i’or  ail  K  .  Then 

lim  An(m)/An*(m)  =  yL(n)/yL(n*)  (4.13) 

m 

for  all  n  i  n*  and  some  n*  ,  0  £  n*  £  o-l. 

If,  In  addition  to  (4.12)  we  are  given  the  series 


1  a  21  Ifcyi(k)  ,  (4.14) 

v.=o 

where 

lim  kKLkyl(k)  -  0  (4.15) 

k— ►» 


for  all  K  ,  then  the  computation  of  y^(n)  by  backward  recursion  based  on 
(3.1)  and  (4.14)  converges. 


Proof: 

Rote  that  y^  ,  if  it  exists,  is  unique,  apart  from  a  constant  multiple, 
by  A- VII.  Since 
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a 


(4.16) 


yh(n) 


nf*  (l  * 

n 


we  have,  by  (3.97) 

\(a)/Tl(»)  a  (ViUJ/y^®))  o(®  h)  >  2  s  h  s  0  > 

far  some  real  crh  ,  and  the  fir*t  part  of  (2.12)  follows. 

Now,  from  (4.15),  we  deduce 

jl^l  a;  Me  **  Ql^kP  ,  k  >  0  ,  0  -  -<*h-2-Re  9L  , 


kwh’vg 

end  we  may  apply  li)  of  Theorem  3.5,  since  (4.12)  implies  that  H0j; 
The  result  is  that  we  can  determine  constants  %  so  that, for  m  > 

K(m;j  sKbyh(m)m’<yh“l/y1(m)  . 

Thus,  by  (4.17) 

lim  RhSh  -  0  ,  2  s  h  *  a 
m— **• 

and  the  proof  is  complete. 


Be  (yiO^OO)  ^Re  *)*!>  , 


(4.17) 

(4.18) 

(4.19) 

.  <  «\>,h  * 
m  , 

(4.20) 

(4.21) 
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Tire  condition  (4-12)  is  rather  stringent,  and  for  a  *  2  or  3  it  can  be* 
weakened  considerably.  We  have,  in  fact. 

Theorem  4.3 

Let  o  =  2  or  3  ,  and  let  jyj00,y2(n) }  be  a  fundamental  set  for  (3.1) 
if  o*2  and  let  | y1(u),y2(n),y3(n)|  be  a  canonical  set  for  (3.1)  if 
o  =  3. 


lim  nL^(ln  n)a  2y^(n^yh(n)  *0,2shic  ,  (4*22) 


n — 


where  L(o)  is  a  positive  number  depending  only  on  a  ,  L(2)  *  0  ,  L(3)  »  1 
Then  for  same  n*  ,  0«n*s  a-1  , 


lim  ^(nO/A^On)  *  y1(n)/y1(n#)  ,  n  *  n*  .  (4.23) 

m— 

Let,  in  addition  to  the  above,  (4,  14)  hold,  with 

1^  -  0  (e^kV*)  (4.24) 

for  some  Q(k)  ,  or  such  that 


eQ(k)kL(o)+ar+i(ln  ) -0(l) 


(4.25) 


t  1 


u 


M 

•■4 


v  «.  u  araMr-**  WM  f*t0p y~i  > 


Then  the  calculation  of  y^(n)  by  backward  re curs ion  bated  on  (?.l)  and 
(4.14)  converges. 


Proof; 


Ve  have,  from  Theorem  3.5, 

Ajj  +  o(l)  ,  £1*1^00  convergent, 

Sh(m) 


0  ^e^BI\na+1-ln  m  y^(m)^  ,  otherwise 


(4.26) 


It  remains  only  to  show  that 


Rj1(m)  »  0  m°~2y^(m^yh(m)}  ,  2  £  h  £  o  .  (4.27) 

Then  (2.12),  and  hence  the  Theorem,  will  follow  from  (4.22)  and  (4.25). 
For  a  ®*  2  ,  (4.27)  is  obvious,  since 


,h-l 


Th(»)  *  (-)  y5.h(«H),  h  «  1,2 


(4.28) 


let,  then,  o  ■  3  .  It  suffices  to  consider  only  !t>  .  Now  |yh(n)j 
Is  canonical,  so,  by  the  construction  of  Theorem  3.3,  we  can  write 


yh(«) 


“  e 


'‘i.h ,  ^.h  i ... 


I .  ia 

In  u 


(In  n  f 


P^(n)  ■  Q^(n)  ♦  @a  In  n  ♦  ln(in  n)  ♦  In 


(4.29) 


(4.30) 
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i 


Since 


[ln(n+l)j  J  -  [la  a]  J  [l  - 


(4.31) 


we  have 


yh(n+l)/yh(n)  « 


*eah  1  +  0  - r) 

'n(ln  n)  / 


(4.32) 


where  An  is  the  forward  difference  operator,  Milne -Thomson  (i960).  Thus 


I^(n-l)  «  -  yI(n)A(n)/[y2(n)B(n)]  , 


A(n)  *  1  -  e 


1+0 


^n(ln  aft)  ] 


B(n)  -  1  -  > .  o(-±—)} 

\n(ln  n)8  /  J 

. i . .4 (v-hy# U , fvv  +  ivV  +  Q(  \l . 

I  a  alnn  VndnnjSyj 


(4.33) 


(4.34) 


(4.35) 


If  f  ,  then  (4.22)  shows  that  pQl  <  p^  and  so  &>(»)— +0  . 
Bence  assume  uei  »  .  If  Uol>  /  Ugj  ,  then  pol  <  Pqj  and  the  leading 

term  of  both  A(n),3(n)  is  unity,  so  again,  by  (4.22),  Rg(m)— *0  .  Hence 
assume  l*0i  -  .  Then  A(n)  is  bounded,  and  B(n)  is  asymptotically 

smallest  in  the  case  where  Qg  *  Qj  ,  Qg  *  ®3  *  ^en  iL  approaches  aero  as 
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l/(n  In  n)  ,  since  we  cannot  have,  furthermore,  pg  -  pj  .  Thus  (4.27) 
follows,  and  R^  ,  R^S^  ,  h  ■  2,3,  approach  zero.  This  proves  the  theorem. 


We  conjecture  that  Theorem  4.3  is  true  for  all  c  ,  i.e.,  that  there 
exists  an  L(o)  for  (4.22)  which  will  Insure  (4.23)  and  (2.7). 

For  o*2  ,  (4.23)  follows  from  a  result  Gautschi  (1961)  proved  far  a 
second  order  difference  equation  with  arbitrary  coefficients,  and  also  for 
this  equation,  Olvcr  (1964)  has  determined  other  conditions  on  y-^yg  which 
will  guarantee  the  convergence  of  Miller's  algorithm,  based  on  (4.14),  to 
yjU)  . 

We  now  present  several  examples  of  applications  of  the  previous  mate¬ 
rial. 

Let 

y(n)  -  J  e’ta'P^^tndt  ,  n  •  0,1,2,...  ,  (4.36) 

0 

p(t)  -  ^  Syt1,  ,0*2  .  (4.37) 

r*i 

A  single  integration  by  ports  shows  that  y(n)  satisfies  the  difference 


equation 


«  X  L  '  - 

y(n)  -  (n+1)  Y  va^n+v)  ■  0  ,  a  •  1  . 


(4,38) 


From  de  Bruljn  (1961  ,  p.  119)  we  have 


y(*0 


jM 


exp 


1- 


1-i 


P^n  ”  +  ^3°  ”  +  *  *“f  Pa"  f <j(o>n)>  a—>«  * 


a  *l 

**a  *  “  Via°  > 


>  (4.39) 


MS 


a 


o-l 

2 


(H 


lv 


-  a 


0-2 


2-1 


When  a  *  2  ,  y(n)  is  essentially  the  parabolic  cylinder  function 
U(n+£,aj/  V2)  ,  the  theory  and  computation  of  which  are  discussed  by  Miller 
(1964). 

Now,  by  Eirkhoff  (1930,  section  2)  ,  we  knew  there  exists  a  fundamental 
set  for  (4.38),  jyh(n)}  ,  whore  has  as  cn  asymptotic  expansion  the  F.8. 
cn  the  right  of  (4.3?)  with  n  replaced  by  neSntfh-l)  ,  Uh  £o  .  Thie 
gives  o  F,S,S.,  so  there  are  no  more.  FVirthermore,  by  Definition  3.6, 

{  yh^n^}  canocical-*  We  can  identify  y(n)  with  y^(n)  . 

We  have 

|yi(n)/yhn(«0|  a  exP  j^ViK"/0^  ^°rh^^[l  *  0(n  ^a)]j  *  (4*4°) 
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fh(y)  ■  sin(Y-hTr/o)  sin(hn/o)  ,  y  »  arg  a<J_;L  . 

For  |y|  <  n  ,  fh  has  zeros  only  at  the  points  hn/o  and  hn/a 
fjj  is  of  one  sign  for  y  between  these  points.  Since 

fh(0)  *  -  sin2(hrr/. )  <  0  , 


we  have 


fh(y)  <  0  >  hrr/a  -n<  y  <  hrr/o 
or 

fhM  <  0  >  <  Y  <  tt/ct  ,  1  S  h  s  a-1  . 

and  so,  for  these  values  of  arg  ,  condition  (4.12)  of  Theorem  4. 
fulfilled. 

Let 

-  run/,)  ±  Vk  • 

k=0 

Since  the  left-hand  side  is  an  entire  function  of  order  o-l  ,  we 
by  Boas  (1954,  p.  10)  , 

W<  C  le  “k/ff  ,kil  , 


(4.41) 

'TT,  SO 

(4.42) 

(4.43) 

(4.44) 

i  is 

(4.45) 

have, 

(4.46) 
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so 


kK|Lk||y1(k)|  s  MlkK+1/°"^  exp  ■  -  |  (l+ln  o^l+Mgk’1^]  =  o(l)  ,  (4.47) 

for  all  K  . 

Also 


3^ky1(k)  =  ]T  VOO  =  1 


k“0  k=0 


(4.48) 


Thus,  by  Theorem  4.2,  the  computation  of  the  integral  (4-33)  by  backward 
recursion  based  on  (4.38)  and  (4.48)  converges  when  |arg  a^-^  |  <  tt/ot  . 

Our  next  example  is  a  class  of  hypergeometric  functions.  (For  notation, 
see  Erdelyi  et  al,  1953;  v.  I.)  Let 


(-An(g*i)n  U 


Q+l 

TT  <aA 

jal  J  /  n+aL^ • • • > n+aQ+2 

q  Q+2*Q+1 n+b^, . . . ,n+bq,2n+Y+ 

TT 

J-l 


Jx)  ’ 


(4.49) 


where  Q,n  are  non-negative  integers,  p+1  ,  y  ,  ,  bj  are  complex  con¬ 

stants,  (aq+2  ■"  3+l)  >  none  of  which  are  negative  integers  or  zero,  and  \ 
is  a  complex  variable,  finite;  £  0  ,  |arg(l-\)  |  <  it  . 


03 


Then  (Wimp  (1966))  y(n)  satisfies  the  difference  equation 


21  [^  +  Ny/\]  y(n+v)  -  0  , 
v«*0 


where 


"  X  >  1*0  •  H^+2  a  0  j  and 


(")  (2n+y )q+3  JJ  (n+Y+v-aj) 

m  « - Jfi - - - 

^  Q+2 

T(Q+3-v)(2n+Y+v)v+I  jj*  (n+a, ) 


/  v-Q-2,2n+Y+v,n+Y+vn-ai, . .  .jn+Y+v+l-a.  ,  v 
Q+4  Q+3\2n+Y+2v+l,n+Y+v-ap...>n+Y+v-aQ+2  \y  >  (4*; 


for  1  S  v  n  Q+2  ,  and 


q  Q 

(-)  (2n*-Y)q+3  TT  (n+Y+v+l-b,) 

Ny - J=i - 1 

Q+2 

r(Q+2-v)(2".+Y+v+l)v  y  (n+Uj ) 

J=1  J 


f  v-Q-l,2n+Y*v+l,ri+Y+v+S-V , . *  .,n+Y+v+2-bc 


i/  '  ^  I 

A  Q>fiFQ+i^2R+Y+2v+i; 


n+Y+v+l-bj_,  y+v+l-bQ 


I1)  ’ 


(4.52 


for  1  s  v  <  Q+i  . 


Note  that  (4.5l)-(4.52)  ore  terminating  hypergeometric  functions,  and 

r°r  any  valu«  *><“  v  >  they  are  ratic.al  functions  of  n  .  Thus  (4.50)  is  of 
the  form  (3.1)  with  tu  a  1  . 


A  result  of  Wimp  and  Luke  (1952,  the  Corollary  on  p.  7,  with  m  =  0  and 
then  u>  =  0  )  shows  that 


l-Z^  7ft)  • 


k=0 


(4.53) 


Since 


(“Q“2)v  ,  -lx 

My  '  — -  +  0(n  )  ,  n— >»  , 


(4.54) 


^  +  0(n  X)  ,  n->» 


(4.55) 


there  are  owo  linearly  independent  solutions  of  (4. 50)  with  the  behaviour 


yx(n)  ~  (-)n[c1(X)]n.l{n)  , 


n— 


(4.56) 


y2(n)  -  (-)n[caU)]  s2(n)  >  n‘ 


(4.57) 


see  Wimp  (1366),  where 


C]U)  3  [l  -  ^]2/x  ,  C2U)  3  [l  +  ^]2/X  7  (4.58) 


and  we  define 


_ _  ,  ,4  |[ars(X-l)-n] 

Vl-X  ■  |l-X|  e  ,  0<  arg(x-l)  <  2rr 


(4.59) 


I 


It  can  be  verified  that  ^  maps  the  X  plane  cut  along  the  real  axis 
from  1  to  *  ,  (i.e.,  |arg(l-x)|  <  n  )  onto  the  interior  of  the  unit  circle 
( X |  <  1  ,  while  maps  the  cut  plane  onto  |x|  >1  .  Thus  for  the  values 
of  arg  X  considered, 


|Ci(*)|  <  1  >  |C2<*)|  >  1  > 


(4,60) 


and  so  (4.12)  bolds  not  only  for  h  ■  2  ,  but  also  when  yh  is  any  purely 
algebraic- logarithmic  solution  of  (4.50).  We  now  proceed  to  show  that  all 
the  remaining  members  of  any  canonical  set  jyfc(n)}  to  which  y^,y g  belong 
are  purely  algebraic -logarithmic. 

Consider 


%(m)  =  (-v)^  >  0  *  m  s 


(4.61) 


%(m)  -  ]T  ('v)mNv  /  0  s  m  s  Q+l  . 


(4.62) 


Using  (4.51)  in  (4.61),  interchanging  the  order  of  summation  and  evaluat¬ 
ing  a  gFi  of  unit  argument,  we  have 


r^(m) 


HTc 

,  .  TT  (a.-Y-m-x) 

.  w?(gnty)Q^3  (Q+8-m)  -  M  J 
Q+2  x  (2x+Y+m)m+ i 

(n+aj  )r(Q+3-a)  J 


(4.63) 


where  A^+2"m)  is  the  fca-ward  difference  operator,  see  Milne-Thomoon 
(i960).  But  the  quantity  in  brackets  above  may  be  decomposed  into  a  poly¬ 
nomial  of  degree  Q+l-ra  in  x  ,  the  indicated  difference  of  which  is  zero, 
and  a  sum  of  terms 


iii 

y~  dj/(2x+  V+m+,J )  . 


(4.64) 


The  dj  may  be  calculated  by  multiplying  the  bracketed  quantity  in  (4.63)  by 
2x+Y+m+J  and  letting  x — »-(Y+m+,l  )/2  .  Then  substituting  (4.64)  in  (4.63) 
gives 


_  (-f* ^  <->r  r(n+ 

TT  '  .  ^  r! 

2  ^n+a.)  V  2  J 


and  an  order  estimate  for  this  sum  is  easily  obtained  by  using 


r(n+o')/r(n>-g )  a  na”^jl  +  0(n  ^)J  . 


(4.66) 


The  result  is 


•Q+2 

1  jj  (Y-2aj)  1  [l  +  0(n_1)J  . 


^(»)  -  (-)Q+v-^4^ 


(4.67) 


liimilarly,  one  finds  that 


Q 

4(-)Q+Wnm'Q^m)<  (v+2-2bj)  *  [l+0(n'l)j,  0 

Oja}  r  ) 

|o(l)  ,  m  •  <1+1  . 


£  m  s  Q  , 


(4.68) 


Next,  we  have 


,k  .  f  i^-4—  4'n) 


(4.69) 


the  B's  being  Bernoulli  numbers,  see  Norlund  (1954,  p.  150). 

We  may  use  the  above  results  to  finally  arrive  at  the  estimate  for 
in  Theorem  3.3, 


4X_1(-)VQA^)  (y+2-2bj )  •  [l+0(n_1)], 


OsksQ, 


-k_ 

n  \  *  .O.]. 

0(n  )  ,  Q+l  s  k  s  Q+2  . 


(4.70) 


»,  in  Theorem  3.3  we  have  p«l,T*-<l,kj*j,k^*Q,  and 


^  ft  <***-*> 

J=0  J  r=0  s=l 


4X'1(-)Q  JT  (Y+0+2-2bs) 

3*1 


(4.71) 


and  oo 


9h  r  2hh  "  Y  -  2  ,  u  =  1,2, . .  .,Q  .  (4.72) 

According  to  that  theorem,  there  are  exactly  Q  algebraic -logarithmic 
solutions  of  (4.50),  and,  if  none  of  the  0^  differ  by  integers,  these  are 
of  the  form 

0  r 

yh+2(n)  -  a  h[1  +  0(n-1)]  ,  lihsQ  .  (4.73) 

Otherwise,  logarithmic  terms  may  appear. 

We  have  thus  determined  a  canonical  set  far  (4.50),  { y^njl  .  A  result 
of  Luke  (1968),  which  is  a  generalization  of  Watson's  result  for  a 2f]_  (see 
Erdelyi  et  al  (1953);  v.  X,  p.  77,  formula  (16)),  enables  us  to  identify 
y(n)  with  a  constant  multiple  of  y^(n)  .  Lastly,  by  examining  as 

given  by  (4.51),  we  see  that  if 

Sj  t  --y-Q-2,  -Y-Q-3, . . . ,  1  S  J  s  QtS  ,  (4.74) 

then  }  0  for  n  s  0  .  If  this  condition  is  satisfied,  as  well  as  the 

conditions  immediately  following  (4,49),  then  Theorem  4.2  may  be  invoked; 


In  particular,  ve  have  demonstrated  a  way  of  computing  Gauss*  function 
wherever  it  is  analytic,  i.e.,  far  |arg(l-X)|  <  it  ,  as  is  seen  by  letting 
Q«0,a^»a,S2*  0+1  *  b,  y+1  ■  c  in  (4.49)-{4.52): 
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s 


9 


(4.75) 


y(n) 


t-ft’H Mn  ,  /n+a,n+b 

(C-l)2n  2  2n+C 


10 


M. 


M, 


(2n+c-l)[2ng+2n(c+l)+c(a+b+l  )-2ab] 
(2n+c)(n+a)(n+b) 


*2 


(n+c+l-a)(n+c+l-b)(2n+c-l)(2n+c)  _  (2n+c-l)(2n+c) 

(n+a)(n+b)(2n+c+2)(2n+c+3)  *  1  (n+a)(n+b) 


(  (4.76) 


To  further  illustrate  the  power  of  the  method,  we  compute  the  function  at  a 
point  on  the  circle  of  convergence  of  its  Taylor  series 

X  =  e"1/3  =*  (1  +  V3i)/2  ,  (4.77) 

with  a  *  2/3  ,  b  »  1  ,  c  ■  4/3  .  It  is  known,  e.g.,  Erdelyi  et  al  (1953, 
v.  I,  p.  106  (55)) 

y(0)-2Fl^2^1|eTTi/3'')=  ■  2nr(l/?i  ^  0.88331  9376  +  0.50998  46791  .(4.78) 

V 4/3 1  y  9r(2/?  f 


Computation  of  y(0)  by  backward  recurrence  using  (4.50),  (4.53)  yields  the 
following  table: 
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a 


ro(«) 


0.80239  8541  +  0.50945  30361 
0.88331  4192  +  0.50998  1S871 


0.88331  9347  +  0.50998  46631 


em  8  y(o)-r0(m) 

(9.2  +  5.31)  x  IQ'4 
(5.2  +  3.01)  x  10“< 
(3.0  +  1.81)  x  10_< 


and  ro(l5)  agrees  with  y(0)  to  all  the  places  (4.78). 

The  series  for  y(0)  converges  only  conditionally*  Knopp  (1947, p.  401)* 
but  this  is  irrelevant*  because  the  Miller  algorithm  will  work  whether  the 
Taylor  series  converges  or  not*  as  long  as  |arg(l-X)  |  <  tt  *  and,  as  is  easily 


seen  from  the  formulation  (2.10)*  the  convergence  is  exponential. 


y(n)  *  rn(m)[l  +  0(®CCl(*)m)] 


,  m— 


(4.79) 


In  this  case 


Cl(a)  -  i(2  ,  |CL(X)|  =  0.260 


•  •  •  • 


(4.80) 


As  our  final  example*  we  take  the  confluent  hypergeometric  function 


yjfa) 


y(n+6*6+i-b*X)  • 


(4.01) 


(Our  notation  and  subsequent  analysis  draw  heavily  on  the  material  contained 
in  Erdelyi  et  al.  (1953,  v.  I*  Ch.  6).) 

y1(n)  satisfies  the  difference  equation 


,  ,  (n+l)[(2n+6+b+l)+x]  ^  fn+l)fii42)  , 

yM  '  I^KiSJj - y(n+l)  +  (5»)U)  y(n+£)  '  0  ’ 


and 


,XM  -  rto.twVMA.-fcW^  ,  0(n-i)j  , 


Jarg  X|  s  it  ,  V(X)  -«fr Tx^6+b^2‘V/2/  f(6)r(b)  , 


1  *  21  »  l“*  x|  <  n  . 

k»0 


(4.82) 

(4.03) 

(4.84) 

(4.85) 


Another  function  satisfying  (4.82)  is 

y2(“)  “  “7^  ?(n+fi,«+l-b;X)  , 

-  r(b)r(6+l-b)V(x)n^6+b^2  "  5/4e2n^X^[l+0(n"^)  j  , 


|  <vg  \  |  s  n 


(4.86) 


(4.87) 


It  is  thus  seen  from  Theorem  4.2  that  the  computation  of  by  backward 
recursion  based  on  (4.8°),  (4.85)  converges  for  jarg  X  |  <  n  .  To  illustrate 
this,  let  6  ■  b  *  £  ,  so 

yx(o)  -  x*Y(ii;x)  -  JW  ,  (4.88) 
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yi(l)  -  ^  Y(3/2,1;X)  -  X  |^X**<l/2,ljX)] 


(4.89) 


«  Jr  yxV/\<x/2)]  -  ^  eX/8[-Kl(A/2)+K0(x/2)(l+l/x)]  •  (4.90) 


For  X  *  4  ,  standard  tables  give 


yx(0)  »  0.94960  6012  ,  yx(l)  «  0.011712616 


(4.91) 


Taking  m  »  10  in  (2.2)  to  (2.6)  yields 


To(l0)  «  0.94961  150B  •••  ,  ^(10)  »  0.01171  2759  ...  ,  (4.92) 


with  approximate  absolute  errors  of  3.3  x  10" 6  and  1.4  x  10"7  ,  respec¬ 
tively.  It  is  interesting  that  the  difference  equation  (4.82)  serves  to 
compute  Ky  ,  while  the  usual  recursion  relation  does  not. 


If  arg  X  *»  tt  ,  then  we  can  define 


<pji(«) 


*£(«) 


»  »e 


C-2i»^|  A 


[l  +  0(m"^)]  # 


■  me 


(4.93) 


(4.94) 


and  Theorems  4.1  and  2.4  apply,  vith  P  -  k  -  2  ,  pg  -  -4jXl  ,  ■  0  ,  and 


(4.95) 


Thus  both  the  functions  (4.81)  and  (4. OS)  can  be  computed  in  this  case, 
provided  that  suitable  normalization  relationships  (or  initial  values)  are 
known,  since  it  is  clear  the  series  (4.85)  will,  in  general,  no  longer  suffice. 

An  even  more  efficient  algorithm  for  the  calculation  of  the  ¥  function 
can  be  based  on  the  third  order  difference  equation  satisfied  by  the  functions 


z(n) 


a 


MgJant*rl.)r(c4a.-l)  31/  jl-n,n+a\ 

r(a)r(6)r(n+b)  ^23  V|  a,6,b  J 


(4.96) 


•j—y  (2n+a-l)r(n+a-l)Xn+a  ~  |\6“aY(n+6,6+l-bj\)  J  ,  (4.97) 


(see  Wimp  (1966,  1967)  and  Luke  and  Wimp  (1963)).  This  is  because  (4.96) 
behaves  as 

d  -3Xl/3n2/3r  -1/3  l 

z(n)  -  c  n  e  ^1  +  0(n  '  )j  ,  n— *■»  ,  (4.98) 

for  some  c  ,  d  (independent  of  n  ),  |arg  \|  <  3n/2  .  (A  canonical  set  for 

2^4 

the  equation  is  readily  obtained  by  replacing  n  above  by  ne  , 
h  -  1,2  . ) 

In  fact,  the  hypergeometric  functions  discussed  in  the  above  three  ref¬ 
erences,  which  are  of  the  farm 

P+11/'  |1-^n+Y+i>bQN 

^Q+2^P+l\Map,g+l  J  •  n  ”  ,  (4.99) 
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can  all  be  computed  by  Miller's  algorithm,  using  the  recursion  relations  and 
normalization  series  given  by  Wimp  (1966),  provided  |arg  X |  <  tt  for 
P  >  Q+i  ,  and  |arg(l+X"i)  j  <  tt  for  P  ■  Q+l  ,  in  which  case  (4.99)  is 
related  to  (4.49).  For  P  >  Q+l  ,  the  difference  equation  for  (4.99)  has  Q 
algebraic- logarithmic  solutions  with  the  same  values  of  9^  as  given  by 
(4.72),  plus  an  additional  P+l-Q  anormal  F.S.E, 

-(m-O^Xe21”*)3^**1"^  ,  h  -  0,1,...,  P-Q  ,(4.100) 

while  if  P  <  Q  ,  then  p  =  oi  =*  1  in  all  the  F.S.S.  of  the  difference  equa¬ 
tion.  Thus  the  recursion  relation  for  (4.99)  has  canonical  sets  whose  com¬ 
ponents  exhibit  widely  varying  behaviour,  depending  on  the  relation  between 
P  and  Q  ,  and  the  equation  may  be  expected  to  furnish  a  number  of  addi¬ 
tional  interesting  applications  of  the  theory  developed  in  this  work. 


yh(n)  ~  chndh  exp 


'  $ 

**■**:•  -^u'y  "  •■".% r*fvr*t,  '****$*? .-‘.'f+xu't * 


APPENDIX 

Here  ve  set  forth  the  notation  used  in  the  body  of  the  thesis,  and  list 
some  frequently  invoked  results  from  the  theory  of  linear  difference  equa¬ 
tions. 

For  determinants  and  matrices,  ve  use  the  notation 


“11 

“12  *  *  * 

“IT 

“21 

• 

• 

“22  ‘  ' 

• 

• 

“2T 

• 

• 

“t1 

• 

“t2  ‘  ‘  ' 

• 

“tt 

(A'l) 


v]I 


“11  “12 


“21  “22 


«Tl  “t2  '  ’ 


(A.2) 


a.  .  always  denoting  the  element  in  the  jth  row  and  hth  column  of  the  de- 
j,n 

terminant  or  matrix. 

Let  <pn  ,  be  (complex -valued)  functions  of  n  »  0,1.2,...  .  We  write 


?n  "  CW  *  n“5'“  » 

if  $n  j*  0  far  n  >  nQ  and  if  a  positive  A  exists  for  vh.ch 


(A.3) 
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If  AH  A  ,  ■>  >  nc 


and 


<Pn  “  o(^n)  ,  n-*«  , 


if,  given  c  >  0  ,  there  exists  an  N(e)  2  nQ  such  that 


<pn/i|(nj<  c  for  all  n  >  N 


Let  p  be  an  integer  2  1  . 


9n  ~  ^n[co+cln”l/P+c2n’2/P+  '  »  c0  f*  0  ,  n 


means  that 


(A.4) 


(A.w') 


(A.6) 


(A,  7) 


V*n  ’  X  Chn_k/P  I  =  <^(r+l)/P)  '  n”>w  ' 

k=0  1 


(A.8) 


for  each  r  =  0,1,2,... 

Occasionally  9  ,  ijt  will  be  functions  of  a  complex  variable  X  belong¬ 
ing  to  some  sector  S  in  the  complex  plane.  Then 


<pU)  *  o(j(xj),  }x|->«  ,  Xes  , 


(A.  9) 


etc.,  ore  interpreted  similarly.  See  Erdelyi  (1956)  or  de  Bruijn  (1961)  for 
details. 
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(A.  13) 


y(k+v)  »  orv  ,  v  =  0,l,2,...,o-l  . 

A-IV.  Equation  (A.  12)  possesses  a  linearly  Independent  set  of  solutions 
|  yh(n)| ,  1  s  h  &  a  ,  called  a  fundamental  set,  and  any  solution  of  (A. 12), 
such  as  (A.  13),  can  be  expressed  as  a  linear  combination  of  these  functions. 

A-V.  Let 

1° 

D(n)  a  |  yh(n+J-l)|  •  (A. It) 

Then 

D(n+1)  -  (-)°D(n)/C0(n)  ,  n  *  0,1,2,...  .  (A. 15) 

A- Vi.  The  equation 

q_ 

/>  C0_v(n+v)y*(n+v)  *  0  ,  n  =  0,1,2,...  ,  (A.lfi) 

v=0 

is  called  the  equation  adjoint  to  (A. 12),  and  the  functions  Th(n)/D(n)  , 

1  s  h  s  o  are  linearly  independent  and  satisfy  (A. 16)  where 

yL(n+l), •••  yh.1(n+l),yh+l(n+l), •• -,ya(n+l) 

e 

\(n)  »  (-)h“l  *  .  (A.lv; 

e 

y^n+o-l),  •  •  •  y^U+o-l^y^n+o-l),  •  •  ->ya(n+a-l) 


Thus 


o 


Th(n)yh(n+r)/D(n) 


l,  r  *  0  , 

0,  1  s  r  £  a-1  , 
-l/C0(n)  ,  r  »  a  . 


A-VII.  Let  the  functions 


{yh(n)}  > 


s  h  s  o  ,  be  such  that 


lim  yh(n)/yh+l(n)  a  0  ,  1  s  h  s  o-l  . 
n— 


Then  the  functions  are  linearly  independent. 


(A.  10) 


(A. 19) 
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